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Abstract 

Based on 58 million J/ip events collected with the BESII detector at the BEPC, the baryon 
pair processes J/tp — > and J/ip — > are observed for the first time. The branching 

fractions are measured to be B(J/ip -> E+E - ) = (1.50 ± 0.10 ± 0.22) x 10" 3 and B(J/ip -> 
= (1.20 ± 0.12 ± 0.21) x 10~ 3 , where the first errors are statistical and the second ones 
are systematic. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.20.Gk, 14.40.Cs 
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I. INTRODUCTION 



Since the discovery of the charmonium states J/ip and ip(2S), a number of baryonic 
decay d>an„e. S have been stud.ed by severa. different experiments |1 fl flfl fl QM 
|9|, [10J, [III, [12j. Baryon-antibaryon decays, which provide a test of the predictive power 
of QCD, have especially attracted interest of both theoretical and experimental experts. 
Among these decays, J/ip — *> and J/ip — > have been studied by DM2 and 

Markll (3, O], but their isospin partners decays J/ip — > and J/ip — > have not 

been measured before. In this article, we study J/ip — > and J/ip — > using the 

large J/ip data sample accumulated with the BESII detector. The decay mode J/ip — > 
is not studied here since the final states contain a neutron and an antineutron, which 
are difficult to detect with the BESII detector. Isospin invariance predicts £>(S + S - ) = 
and £>(H°H°) = However in the quark model, there are well-known 

nn 

isospin breaking contributions in J/ip baryonic decays [131. Il4j|. 



II. THE BESII DETECTOR AND MONTE CARLO SIMULATION 

BESII is a conventional solenoidal magnet detector that is described in detail in 
Ref. fl5| . A 12-layer vertex chamber (VTC) surrounding the beam pipe provides trigger 
and track information. A forty- layer main drift chamber (MDC), located radially outside 
the VTC, provides trajectory and energy loss (dE/dx) information for charged tracks over 
85% of the total solid angle. The momentum resolution is <J p /p = 0.017-\/l + p 2 (p in 
GeV/c), and the dE/dx resolution for hadron tracks is ~ 8%. An array of 48 scintillation 
counters surrounding the MDC measures the time-of-flight (TOF) of charged tracks with 
a resolution of ~ 200 ps for hadrons. Radially outside the TOF system is a 12 radiation 
length, lead-gas barrel shower counter (BSC). This measures the energies of electrons and 
photons over ~ 80% of the total solid angle with an energy resolution of cte/E = 22%/y/E 
(E in GeV). Outside of the solenoidal coil, which provides a 0.4 Tesla magnetic field over 
the tracking volume, is an iron flux return that is instrumented with three double layers 
of counters that identify muons of momentum greater than 0.5 GeV/c. 

n 

In this analysis, a GEANT3 based Monte Carlo (MC) simulation program [16j with 
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detailed consideration of real detector responses (such as dead electronic channels) is 
used. The consistency between data and MC simulation has been carefully checked in 
many high-purity physics channels, and the agreement is quite reasonable [171 ]. 



III. EVENT SELECTION 

The data sample used for this analysis consists of (58.0±2.7) x 10 6 J /if} events collected 
with the BESII detector [3]. The decay channels investigated are J /if) decays into 
and baryon pairs, where S + decays to ir°p (tt° — > 77), S° to tt°A (A — > n^p). 

Therefore the final states for the two decays are ppyyyy and 7r + 7r~pp7777, respectively. 
Both decays contain four photons in the final states. Candidate events are required to 
satisfy the following common selection criteria: 

1. Events must have two or four good charged tracks with zero net charge. A good 
charged track is a track that is well fitted to a helix in the MDC, and has a polar 
angle, 9, in the range |cos#| < 0.8. For J /if) — > E + E~, tracks are required to 
originate from the interaction region of R xy < 0.02 m and \z\ < 0.2 m, where R xy 
is the distance from the beamline to the point of closest approach of the track to 
the beamline, and |z| is the distance along the beamline to this point from the 
interaction point. For J/if> — > because of the long lifetime of 5 and A, tracks 
are not required to originate from the interaction region. 

2. The TOF and dE/dx measurements of the charged tracks are used to calculate 
Xpid va hies for the hypotheses that the particle is a pion, kaon, or proton. Only 
the two proton tracks must be identified with the requirement that Xpid f° r the 
proton hypothesis is less than those for the 7r or K hypotheses. 

3. Isolated photons are those that have energy deposited in the BSC greater than 50 
MeV, and the angle between the photon entering the BSC and the shower devel- 
opment direction in the BSC is less than 37°. In order to remove the fake photons 
produced by p annihilation and those produced by hadronic interactions of tracks, 
the angle between the photon and antiproton is required to be larger than 25° and 
those between the photon and other charged tracks larger than 8°. 




FIG. 1: Invariant mass distributions of two photons for J/ijj — > candidate events, (a) 

Scatter plot of M(jj)i versus -^(77)2 for the combination with the minimum R(jr°), (b) dis- 
tribution of M(77)i, (c) distribution of -^(77)2 in the data sample, and (d), (e), and (f) are 
the corresponding plots for MC simulated J/ip events. 

Different kinematic fits are used in the selection of the two decay channels. A four 
constraint (4C) kinematic fit under the PP7777 hypothesis is performed for J/ip — > 
If there are more than four photon candidates in an event, all combinations are tried, and 
the combination with the smallest x\c * s retained. We require the minimum x\c t° t> e 
less than 15. For J/ip — > a six constraint (6C) kinematic fit under the hypothesis 

J /if) — > 77777r + 7r _ pp with the invariant mass of the two photon pairs constrained to the 
7T° mass is performed, and the x 2 °f the 6C fit is required to be less than 50. 
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IV. DATA ANALYSIS 

A. J /if) -» E+E - 

The candidate events for this decay mode contain two 7r°, and there are three pos- 
sible combinations of (77)1 (77)2 to form a ir° pair. The 7r° pair with the minimum 
R(n°), where R(ir°) = ^(M (77)1 - M(vr )) 2 + (M( 77 ) 2 - M(vr )) 2 , is chosen for further 
analysis. Figure Q] shows the mass distributions of candidate events with the minimum 
i?(7r°) for data and MC samples, respectively. A clear 7r°7r° signal is observed in the 
data sample. In order to select tt° pair events, |M(77)i — M(tt )\ < 0.03 GeV/c 2 and 
|M(7 7 ) 2 - M(vr )| < 0.03 GeV/c 2 are required. 

After 7T° selection, there are still two possible ir p combinations from which to form 
the £ of the £ + £~ pair. The combination having the smallest value of i?(E) = 



J(M(tt° {1) p) - M(£)) 2 + (M(n° {2) p) - M(£)) 2 is selected for further analysis. Figure M 
shows the up invariant mass for data and MC simulated J/ip —>■ events. A clear 

signal is seen in the bottom left corner of Fig. |2^a). Figure EJ^b) is the M(ir^p) 
distribution by requiring \M{n° {2) p) - M(E)| < 0.03 GeV/c 2 , (c) is the M(ir° {2) p) distribu- 
tion by requiring |M(7r^p) — M(E)| < 0.03 GeV/c 2 , and (d) is the sum of (b) and (c) 
scaled by a factor of 0.5. 

Possible backgrounds come from channels with pp production, including J/ip — -> pp, 
J /ip —>■ jpp, J/ip — * pp7t°, J ftp — > ppi] (r] — > 77 or 7] —>■ 3tt°), J/ip —>■ ppu (u — > 77T°), and 
J /ip — > pp7c°TT . MC events are generated with a phase space generator for the first two 
background channels. None or only a few events survive the selection criteria; therefore 
contamination from the first two channels is negligible. For J/ip —>■ ppir , J /ip — > ppi], and 
J/ip —> ppu, MC events are also generated according to phase space. Using the branching 
fraCtion8 fr0m the PDG Q, the nUmberS of events fr0m the 8 e ehaoneis are expected to 
be 7.2, 72.4 and 7.8 in the whole np mass region, respectively. For J '/ip — > ppn ^, the 
branching fraction is unavailable, so the normalized number of events for this background 
can not be determined. However, the np invariant mass distribution from all backgrounds 
is smooth, so these backgrounds will not affect the determination of the number of signal 
events in fitting the irp mass distribution. 
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FIG. 2: Distributions for J/ij) — » S + S candidate events, (a) Scatter plot of M(tv^p) ver- 
sus M(ir° (2) p); (b) M(i$wp) distribution by requiring \M(ir%,p) - M(S+)| < 0.03 GeV/c 2 , (c) 
M(7r| ) 2 - ) p) distribution by requiring \M (n^p) — M(S + )| < 0.03 GeV/c 2 , and (d) is the sum of 

(b) and (c) scaled by a factor of 0.5. (e), (f), (g) and (h) are the corresponding plots of (a), (b), 

(c) and (d) for MC simulated J/tp — ► events. 

In order to determine the branching fraction, we fit the S signal in Fig. [21(d) with a 
histogram of the signal shape from MC simulation together with a second order polynomial 
for the background. The fit is shown in Fig. [31 and it yields 399 ± 26 signal events, with 
the goodness of the fit being x 2 /ndf = 22.4/31 ~ 0.72. The J/ip — > efficiency is 

determined to be e—1.75% using MC simulated signal events, and the branching fraction 
is, 

B{ JH, -+ £+2-) - N <®'< 



N(J/ip) ■ £ 2 (£+ -> 7T°p) ■ B 2 (ir° -> 77) 
= (1.50 ±0.10) x 10^ 3 , 

where the error is statistical. 
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FIG. 3: Fit to np invariant mass of J/ij) — > candidate events with MC simulated signal 

shape and a second order polynomial as background shape. The shaded histogram is background 
from MC simulated J/ip — ► ppn , J/ifi — > ppr], and J/tp — > ppuj, normalized according to the 
branching fractions in PDG [19|]. The dashed histogram shows the shape of MC simulated 
J/ip — > ppir°Tr normalized using an assigned branching fraction of 0.003. 

B. J/%b -» 

The candidate events for this decay mode contain a AA pair. In order to select 
AA events, we require the n~p and 7i + p invariant masses satisfy \M(ir~p) — M(A)| < 
0.01GeV/c 2 and \M(n + p) - M(A)| < 0.01GeV/c 2 . 

There are three possible combinations of the four photons to form a tt° pair in 
the 6C kinematic fit; the combination with the minimum xlc ^ s considered to be the 
correct one and selected for further investigation. After 7r° pairs are selected, there 
are two possible combinations (tt^A, 7rfLA) to form 5° candidates. Analogous to 
the analysis of S + E~, we choose the combination with the lowest value of i?(H°) = 
^(M(tt° (1) A) - M(H )) 2 + (M(7rj ) 2) A) - M(S )) 2 for further study. Figure 1 shows the 
7rA invariant mass for this case for data and MC simulated J/ip — > events. In Fig. H] 
(a), besides the clear S(1385)°S(1385)° production is also visible. In Fig. H] (b) is 
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the M(?r? 1) A) distribution after requiring |M(7rLA) - M(H°)| < 0.03 GeV/c 2 , (c) is the 
M(tt[ ) 2) A) distribution requiring \M(n° {1) A) - M(E°)\ < 0.03 GeV/c 2 , and (d) is the sum 
of (b) and (c) scaled by a factor of 0.5. 

Possible backgrounds come from channels with A or S production, including J /if) — > 
J/ip -> S°7r A(+c.c), J/^ -»■ £(1385)°£(1385_) , and J/^ -> 7r°7r°AA. Using the 



branching fraction of J/^ — > S°S° from the PDG 19] and assuming isospin invariance 
holds for J/ip -> SttA + c.c. and J/V> -> E(1385)S(1385), we obtain 0.6, 45.9, and 51.3 
background events from MC simulation, respectively. The shaded part in Fig. [5] shows 
the normalized simulated backgrounds from the first three channels, which do not exhibit 
any peaking structures in the S° mass region. Since the branching fraction is unavailable 
for J /if) — > 7r°7r° AA, we do not determine the normalized number of events for this de- 
cay mode. However, MC simulation indicates that the 7rA invariant mass distribution is 




(b) s 



(c) ~- 



M(ti A) (GeV/<0 




(g) c 



M(itA) (GeV/O 



(h) 



M(ji"a) (GeV/c*) 



M(ji"a) (GeV/cT 



M(jt"A) (GeV/c^) 



M(iA) (GeV/c^) 



FIG. 4: Plots for J/ip — ► candidate events, (a) is the scatter plot of M(7r? 1 \A) versus 

M(7t9 2 n A), (b) is the M(7nLA) distribution recoiling against 3°, selected by requiring \M(n^A) — 
M(H°)| < 0.03GeV/c 2 , (c) is the M(7T/L A) distribution recoiling against H , selected by requiring 
| M A) — Af (H°) | < 0.03GeV/c 2 , (d) is the sum of (b) and (c) scaled by a factor of 0.5, (e), (f), 
(g) and (h) are the corresponding plots of (a), (b), (c), and (d) for MC simulated J /if) — > 
events. 
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smooth without any peaking, and therefore this background will not affect the determi- 
nation of the number of signal events. We also studied backgrounds from AA sidebands 
and other possible background channels listed in the PDG, but the contaminations were 
found to be negligible. The fitted number of signal events is insensitive to the shape of 
all the backgrounds considered here. The numbers of signal events using different back- 
ground shapes differ slightly from each other; we consider this difference as one source of 
systematic error. 

To get the number of signal events, we fit the observed 5 signal in Fig. H] (d) by a 
histogram of the signal shape from MC simulation plus a second order polynomial as the 
background. The fitting result is shown in Fig. and the fit yields 206 ±20. The efficiency 
is determined to be £=0.74% using MC simulated signal events generated according to a 
phase space distribution. The branching fraction is, 

N(E)/e 



B(J/1> 



N{J/lp) • B 2 (E° -> 7T°A) ■ B 2 (A -> TT-p) ■ £ 2 (7T° 

fl.20±0.12) x 10~ 3 , 



77) 



o 

> 




LU 
> 
LU 




M(;tA) (GeV/cr) 



FIG. 5: Fit to irA invariant mass spectrum of J/tp — ► candidate events. Dots with 

error bars are data, the hatched histogram is the normalized background from all the channels 
considered in the text, and the solid histogram is the fit to data using a histogram of the signal 
shape from MC simulation plus a second order polynomial for background. 
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where the error is statistical. 



V. SYSTEMATIC ERRORS 



The systematic errors on the branching ratios mainly arise from the uncertainties in the 
MDC tracking, particle identification, photon efficiency, angular distribution parameter 
a in event generators, kinematic fitting, background shapes, and the total number of J / if? 
events. The errors from different sources are listed in Table [H 

The uncertainties caused by MDC tracking and particle identification (PID) are esti- 
mated by the difference of the selection efficiency of proton and antiproton between data 



and MC simulation 



201 ] . The efficiencies of PID and track reconstruction for protons and 



antiprotons that enter the detector being reconstructed and identified are measured using 
samples of J /if? — ► ir + ir~pp, which are selected using PID for three tracks, allowing one 
proton or antiproton at a time to be missing in the fit Q. It is found that the efficiency 
difference of one proton identification is about 1% for J/if? —>■ and about 2% for 

J /if) — > depending on the momentum of the final state particles. The 7r ± tracking 
and PID efficiencies are simulated within 1% per track. Therefore we get a total of 2% 
for J /if? — > and 6% for J /if? — > respectively. 



TABLE I: Summary of systematic errors 



Source 




J/i/j E°E° 



MDC tracking and PID 
Photon efficiency 
Kinematic fit 
Background shape 
J /if) statistics 
a 

VC trigger efficiency 
Total 



6 

8 

8.4 
5 

4.7 

7 

4 
16.9 
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The photon detection efficiency is studied using J /if) — * p°7r° in Ref. [21] . The results 
indicate that the systematic error is about 2% for each photon. Therefore, 8% is taken 
as the systematic error of photon efficiency for the two decay modes. 

The angular distribution of the baryon in J /if) decay is 1 + a cos 2 6, with 9 being the 
polar angle of the baryon in J/if> rest frame. To estimate the uncertainty originating from 
the angular distribution parameter a, we generate MC samples for a = and a = 1, 
separately. The differences of efficiency between a = and a = 1 are taken as systematic 
errors, which is 8% for J /if) — > E + E~ and 7% for J /if) — > 

The systematic errors for the kinematic fits are 4% for and 8.4% for they 

are taken from earlier studies 22|, |23j. The VTC trigger efficiency systematic errors are 



estimated to be 1% for S + S and 4% for 

The systematic error of the background shape used is estimated by measuring the 
difference of the numbers of fitted signal events for different background shapes. For 
J /if) — > we also fit the irp invariant mass distribution using the normalized back- 

ground, or MC simulated J/if> — ► ir 7r°pp events as the background shape. The difference 
of the numbers of signal events is about 6% compared to our nominal fit, which is taken 
as the systematic error of the background shape. For J /if) — > we fit the 7rA invari- 

ant mass distribution using the normalized background or MC simulated J /if) — > 7r°7r°AA 
events as the background shape, and we estimate a systematic error of about 5%. 

Uncertainty on the total number of J /if) events is 4.7% 18] . Combining these errors 
in quadrature gives total systematic errors of 14.4% for J /if) — > and 16.9% for 

J/if) -> 



VI. RESULTS AND DISCUSSION 



Based on 58 x 10 6 J /if> events accumulated at BESII, we report first measurements of 
the branching fractions of J /if) decays into the baryon pairs S + E~ and The results 

are listed in Table HT1 including the results of J /if) — > TPTP and J /if) — > We note 

that the isospin partners, S + and S° and also S° and H~, have similar branching fractions 
in agreement with expectations of isospin symmetry. Furthermore, according to the phase 
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space corrected branching fraction |Mj| 2 = B(J/ip — > BiB/) / (np* / y/s) |5j, we obtain the 
phase space corrected branching fractions (1.50 ± 0.10 ± 0.22) and (1.45 ± 0.15 ± 0.26) for 
J/ip — > and J/ip — > respectively. We note that the increase of strangeness 

does not greatly change the branching fraction, indicating the flavor symmetric nature of 
gluons. We also calculate the ratios of the ip(2S) results in Ref. [6J] to those from J/ip 
measurements in this article after removing the phase space factor, and obtain (14.3 ± 
5.2)% for S+S- and (17.3 ± 6.7)% for They agree with the so called "12% rule" 

predicted by perturbative QCD [24| within la. 
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TAB 



E II: Branching fractions (10 3 ) of J /if) decays into E + S , and 



111 ]. The first error is statistical and the second systematic. For the results of MarkI [1 




the statistical and systematic errors were added in quadrature. 



Channels 


J/ip -» 


J/iP -> z-z+ 


J/ip -> 




MarkI 


1.3 ±0.4 


1.4 ±0.5 




3.2±0.8(H°H°±H-H+) 


Markll 


1.58 ±0.16 ±0.25 


1.14 ±0.08 ±0.20 






DM2 


1.06 ±0.04 ±0.23 


1.4 ±0.12 ±0.24 






BESII 


1.33 ±0.04 ±0.11 




1.50 ±0.10 ±0.22 


1.20 ±0.12 ±0.21 


BaBar 


1.15 ±0.24 ±0.03 
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